Evidence for the immobile bipolaron formation in the paramagnetic state of the 

magnetoresistive manganites 
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Recent research suggests that the charge carriers in the paramagnetic state of the magnetoresistive 
manganites are small polarons. Here we report studies of the oxygen-isotope effects on the intrinsic 
resistivity and thermoelectric power in several ferromagnetic manganites. The precise measurements 
of these isotope effects allow us to make a quantitative data analysis. Our results do not support a 
simple small-polaron model, but rather provide compelling evidence for the presence of small immo- 
bile bipolarons, i.e., pairs of small polarons. Since the bipolarons in the manganites are immobile, 
the present result alone appears not to give a positive support to the bipolaronic superconductivity 
theory for the copper-based perovskites. 



The discovery of "colossal" magnetoresistance (CMR) 
in thin films of Rei_ x A x Mn03 (Re = a rare-earth ele- 
ment, and A — a divalent element) jjj has stimulated 
extensive studies of magnetic, structural and transport 
properties of these materials [0. The physics of man- 
ganites has primarily been described by the double- 
exchange (DE) model ||. However, Millis, Littlewood 
and Shraiman [0 pointed out that the carrier-spin in- 
teraction in the DE model is too weak to lead to car- 
rier localization in the paramagnetic (PM) state, and 
thus a second mechanism such as small polaronic ef- 
fects should be involved to explain the observed resis- 
tivity data in doped manganites. Following this original 
idea, more theoretical models were proposed (||| ■ On the 
other hand, Alexandrov and Bratkovsky H have recently 
shown that the essential physics in manganites should in- 
volve the formation of small bipolarons in the PM state 
in order to explain CMR quantitatively. 

The first experimental evidence for small polaronic 
charge carriers in the PM state was provided by trans- 
port measurements It was found that the activa- 
tion energy E p deduced from the conductivity data is 
one order of magnitude larger than the activation en- 
ergy E s obtained from the thermoelectric power data. 
Such a large difference in the activation energies is the 
hallmark of the small-polaron hopping conduction. The 
giant oxygen-isotope shifts of the ferromagnetic transi- 
tion temperature Tq give clear evidence for the presence 
of polaronic charge carriers in this system |j 10 1. More- 
over, the fast and local techniques have directly shown 
that the doped charge carriers are accompanied by local 
Jahn- Teller distortions |il)-|i"4|]. However, all these ex- 
periments cannot make a distinction between small po- 
larons and small bipolarons since both are dressed by 
local lattice distortions. Small bipolarons are normally 



much heavier than small polarons, and should be local- 
ized in the presence of small random potentials. In or- 
der to discriminate between polarons and bipolarons and 
to place constraints on the CMR theories, we studied 
the oxygen-isotope effects on the intrinsic resistivity in 
the high-quality epitaxial thin films of Lao.7sCao.25Mn03 
and Ndo.7Sr .3Mn03. We also measured the thermoelec- 
tric power for the oxygen-isotope exchanged ceramic sam- 
ples of Lao.75Cao.25Mn03. The data cannot be explained 
by a simple small-polaron model, but are in quantitative 
agreement with a model where the formation of small 
immobile bipolarons is essential. 

The epitaxial thin films of Lao.75Cao.25Mn03 (LCMO) 
and Nd .7Sr . 3 MnO3 (NSMO) were grown on <100> 
LaA103 single crystal substrates by pulsed laser depo- 



sition using a KrF excimer laser |15|. The film thickness 
was about 190 nm for NSMO and 150 nm for LCMO. Two 
halves were cut from the same piece of a film for oxygen- 
isotope diffusion. The diffusion for LCMO/NSMO was 
carried out for 10 h at about 940/900 °C and oxygen 
pressure of 1 bar. The 18 0-isotope gas is enriched with 
95% 18 0, which can ensure 95% 18 in the 18 thin films. 
The ceramic 16 and 18 samples of Lao.7sCao.25Mn03 
were the same as those reported in Ref. ff6f . The resistiv- 
ity was measured using the van der Pauw technique, and 
the contacts were made by silver paste. The measure- 
ments were carried out in a Quantum Design measuring 
system. The thermoelectric power was measured using 
an apparatus modeled after a seesaw technique flTj . The 
absolute uncertainty is less than 0.25 /xV/K and the sys- 
tematic error is ±0.1/iV/K. 

Fig. 1 shows the resistivity of the oxygen-isotope ex- 
changed LCMO films above 1.1T C , where T c = 231.5 
K for the 16 sample and 216.5 K for the 18 sample 
(see Table 1). The oxygen-isotope shift of Tc is 15.0(6) 
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K in the films, in excellent agreement with the results 
for the bulk samples Jll|. From the figure, one can see 
that there is a large difference in the intrinsic resistivity 
between the two isotope samples. Such a large isotope ef- 
fect is reversible upon the oxygen isotope back-exchange. 
We should mention that the intrinsic resistivity of the 
compounds can be only obtained in high-quality thin 
films and single crystals. Our LCMO films even have 
lower residual resistivity p than the corresponding sin- 
gle crystals Jl8|] , indicating a high-quality of the films. 
We also checked that the resistivity for a ceramic sam- 
ple of Lao.66Bao.34Mn03 is very different from that for 
the corresponding single crystal; the activation energy 
E p for the former is about 3 times larger than for the 
latter. The temperature dependence of the resistivity 
obtained in the ceramic samples does not represent the 
intrinsic behavior of the bulk, and thus one cannot use 
ceramic samples to study the isotope effect on the in- 
trinsic resistivity. Moreover, the van der Pauw technique 
is particularly good to precisely measure the resistivity 
difference between oxygen-isotope exchanged films whose 
thicknesses are identical. Thus the data shown in Fig. 1 
represent the first precise measurements on the intrinsic 
resistivity of the two isotope samples. 
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FIG. 1. The resistivity of the oxygen-isotope exchanged 
films of Lao.75Cao.25MnC>3. The maximum temperature of 
the data points for the ie O film is 380 K. The solid lines are 
the fitted curves by Eq. 1. As in Ref. ^j, we excluded the 
data points below \.\Tc for the fitting. 

In the PM state, the dominant conduction mechanism 
is thermally activated hopping of adiabatic small po- 
larons. When T>W p /kB (where W p is the polaron band- 
width), the temperature dependence of the resistivity is 
given by [0(2011 



p = BTexp[(E a + E s )/k B T], 



(1) 



where B is a coefficient, depending on a characteristic 
optical phonon frequency u> (i.e., B cx 1/uj )] E s is the 
energy required to excite a polaron from a localized im- 
purity state (see Fig. 2a); E a — (rjE p /2) — t (here E p is 



the polaron binding energy, t is the "bare" hopping in- 
tegral, r\ = 1 for Holstein polarons jL9) and r\ ~ 0.2-0.4 
for Frohlich polarons 21 1). In the harmonic approxima- 
tion, E p and E a are independent of the isotope mass M. 
Moreover, since E s does not depend on M either [fj0| , 
one expects that E p = E a + E s should be independent 
of M. 

The thermoelectric power is given by [p|.p0[ 



S = E s /eT+S 



(2) 



where S is a constant depending on the kinetic energy of 
the polarons and on the polaron density [^0| . One should 
note that Eq. 2 is valid only if there is one type of carriers 
(e.g., holes). Comparing Eq. 1 and Eq. 2, one readily 
finds that E p — E a + E s > > E s . Jaime et al. [|j used the 
above equations to fit their resistivity and thermoelectric 
power data, and found that E p >> E s . This is consistent 
with the small polaron hopping mechanism. 

We now fit our data by Eq. I (see solid lines in Fig. 1). 
It is apparent that the fits are quite good for both isotope 
samples. However, the isotope dependencies of the pa- 
rameters B and E p are not expected from Eq. 1. Upon 
replacing 16 O with 18 O, the parameter B decreases by 
20%, and the parameter E p increases by 11 meV. This is 
in contradiction with Eq. 1 which predicts that the pa- 
rameter B should increase by about 6% while E p should 
not change. Since our measurements are very accurate, 
the unusual isotope effect on B cannot be caused by the 
experimental uncertainty. One possible explanation is 
that Eq. 1 does not hold below 380 K (the maximum 
temperature of the data points in Fig. 1) since the con- 
dition T>W p /kB may not be satisfied in this system. 
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FIG. 2. A schematic diagram of the polaron band, and po- 
laron trapping into impurity (IP) states (a), or into localized 
bipolaron (BP) states (b). 

One should modify Eq. 1 if T<W p /kB- It is known 
that the resistivity can be generally expressed as p — 
l/cr= 1/nejjL, where n is the mobile carrier concentration 
and p is the mobility of the carriers. For adiabatic small- 
polaron hopping, the mobility is given by |ll| 
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^exp(-E a /k B T), 



(3) 



where d is the site to site hopping distance, which is 
equal to a/v2 in manganites since the doped holes 
in this system mainly reside on the oxygen sites [p3| . 
Here E a should also be modified at low temperatures 
§§, i.e., E a = ( V E p /2)f(T) - t, where f(T) = 
[tanh(?iw /4A: B r)]/(?iw /4fc B T) for T > Tiuj /Ak B ~ 200 
K The mobile polaron density n can be easily 

calculated with the help of Fig. 2. When T>W p /k B , 
n ex exp(— E s /k B T) ^0|, which leads to Eq. 1 by com- 
bining with Eq. 3. On the other hand, when T<W p /k B , 
n = 2(2Trm p k B T/h 2 ) 3 / 2 exp(— E s /k B T) if we assume a 
simple parabolic band 0,^O|. Here m p is the effective 
mass of polarons and related to W p by m p — 6h 2 /a 2 W p . 
In fact, the above n(T) expression is the same as that for 
semiconductors when the chemical potential is pinned to 
the impurity levels. Using the above n(T) expression in 
the case of T<W p /k B and Eq. 3, we finally have 
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P = —f=exv{Ep/k B T), 
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(4) 



where C = (a/i/e 2 \/fc^)(1.05iy p ) 1 - 5 /^ o . The quantity 
C should strongly depend on the isotope mass M and de- 




FIG. 3. The resistivity of the oxygen-isotope exchanged 
films of (a) Lao.rsCao.asMnOs; (b) Ndo.ySro.sMnOs. The 
solid lines are fitted curves by Eq. 4. 

crease with increasing M. This is because W p de- 
creases strongly with increasing M according to W p = 
W exp(-TE p /huj ) = IT exp(- ff 2 ) MM. Here W a is 



the bare bandwidth (W — 12i), T is an isotope inde- 
pendent constant, which is less than 0.4 for both Hol- 
stein and Frohlich polarons when the coupling constant 
A = 2E p /W o <0.5 and hu /t<\ ||. Clearly, the pa- 
rameter C should decrease significantly with increasing 
oxygen-isotope mass, in contrast to the parameter B 
(oc 1/uJo) which would increase by about 6% upon re- 
placing 16 O with 18 O. 

Now if one considers that small polarons are bound into 
localized bipolarons (i.e., pairs of small polarons) 0, the 
only modification for both Eq. 2 and Eq. 4 is to replace 
E s by A/2 (see Fig. 2). Here A is the bipolaron binding 
energy and given by A = 2(1 — T)E p — V c — W p , where 
V c is the Coulombic repulsion between bound polarons 
pjj . It is apparent that A and thus E s should depend 
on the isotope mass M due to the fact that W p is a strong 
function of M. 
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FIG. 4. The deviations between the data of the 16 O film of 
Lao.75Cao.25Mn03 and the fitted curves by Eq. 1 and Eq. 4. 

In Fig. 3, we show the resistivity data for the 16 and 
18 films of both LCMO and NSMO compounds, and fit 
them by Eq. 4. The data can be well fitted by Eq. 4 with 
the parameters summarized in Table I. It is striking that 
the fits with Eq. 4 are much better than the fits with 
Eq. 1, as seen more clearly in Fig. 4 where the deviations 
between the data and the fitted curves are plotted for the 
LCMO 16 film. The fit with Eq. 1 has a large systematic 
deviation (>4%), while the fit with Eq. 4 has a much 
smaller deviation (<1%). It is easy to check that the T 
dependence of the prefactor of the exponential function 
in Eq. 1 and Eq. 4 is not important when E p is large, but 
becomes crucial when E p is small. This can naturally 
explain why the resistivity data in the Lai-ajCa^MnOa 
system can be well fitted by Eq. I, but the data of the 
Lao.7Sro.3Mn03 film are not consistent with Eq. 1 p4[ . 
We speculate that the data of Lao.7Sro.3Mn03 in Ref. |24| ] 
should be in better agreement with Eq. 4. 

From Table I, one can see that the parameter C for 
LCMO decreases by 35(5)%, and E p increases by 13.2(3) 
meV. For NSMO, C decreases by 40(7)%, and E p in- 



3 



creases by 14.2(8) meV. The deduced oxygen-isotope ef- 
fects on the parameter C and E p are in qualitative agree- 
ment with Eq. 4. Since this equation is valid independent 
of whether the small polarons are bound into localized 
bipolarons or to the impurity centers (see Fig. 2 and dis- 
cussion above), one can only draw a conclusion that the 
isotope dependence of the electrical transport in the PM 
state of the manganites agrees qualitatively with small- 
polaron hopping conduction. 

Now the question arises: Can the present data tell 
whether the small polarons are bound into localized bipo- 
larons or to the impurity centers? The clarification of this 
issue can place an essential constraint on various CMR 
theories and on bipolaronic superconductivity theory for 
the cuprate superconductors. If the small polarons are 
bound to impurity centers, there will be no isotope effect 
on E s [^0|. Then, the isotope effect on E p only arises 
from the isotope shift of E a because E p = E a + E s . 
Since E a = (rjE p /2)f(T) - W a /12, only the quantity 
f(T) = [ta,nh{hcj /4k B T)}/(hu; /4:kBT) may depend on 
the isotope mass if the temperature is not so high com- 
pared with fiojo/ks- If we take hco = 74 meV, typical 
for oxides |p2l , we have f(T) = 0.855 and the isotope- 
induced change Sf(T) = 0.0134 at T = 300 K, the mid- 
point temperature of the data points for the LCMO films. 
Using the relation: (r)E p /2)f(T) = W /12 + E p - E Sl 
and E s = 13.2 meV for Lao.75Cao.2sMn03 (see below), 
we find rjE p /2 = 0.55 eV. Here we have used W = 4.9 
eV, as estimated from the relation: W — 6h 2 /a 2 nib and 
mi = 0.61m e (||]. The estimated bare bandwidth is typ- 
ical for the oxygen band and in good agreement with the 
electron-energy- loss spectra p3| . Then SE p = SE a — 
{r)E p /2)8f{T) = 7.4 meV, which is about half the value 
observed. This implies that there must be an isotope ef- 
fect on E s , which is only possible if the small polarons 
are bound into localized bipolarons as discussed above. 
The oxygen-isotope shift of E s in the LCMO compound 
is SE S = SE p - SE a = 5.8(3) meV. 

TABLE I. Summary of the fitting and measured parame- 
ters for the 16 and ls O films of Lao.7sCa .25Mn0 3 (LCMO) 
and Nd .7Sr .3MnO 3 (NSMO). The errors of the parameters 
comes from the fitting and from the van der Pauw measure- 
ment. The absolute uncertainty of the thickness of the films 
was not included in the error calculations since it only influ- 
ences the absolute values of the resistivity. 
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LCMO( itt O) 


231.5(3) 


122(2) 


17.3(5) 


72.8(2) 


LCMO( ls O) 


216.5(3) 


141(2) 


12.9(3) 


86.0(1)) 


NSMO( le O) 


204(1) 


248(4) 


23.2(8) 


78.8(4) 


NSMO( ls O) 


186(1) 


289(4) 


16.2(7) 


92.9(4) 



One can also obtain the isotope shift of E s by mea- 
suring the thermoelectric power for two isotope samples 
according to Eq. 2. In Fig. 5, we plot the thermoelectric 
power S as a function of 1/T for the 16 O and 18 O samples 
of Lao.75Cao.25Mn03. Both Tc's and the isotope shift of 
the ceramic samples [fl6f are the same as those in the cor- 
responding thin films. Since the grain-boundary effect on 
S is negligible, the thermoelectric power obtained in ce- 
ramic samples should be intrinsic. From the slops of the 
straight lines in Fig. 5, we find E s — 13.2 meV for the 
16 sample and 18.7 meV for the 18 0. The isotope shift is 
SE S = 5.5(5) meV, in remarkably good agreement with 
that (5.8 meV) deduced independently from the above 
resistivity data. 
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FIG. 5. The thermoelectric power S(T) of the oxy- 
gen-isotope exchanged ceramic samples of Lao.7sCao.25Mn03. 

We can use the values of the parameter C for the 16 O 
films to calculate the polaron bandwidth W p according 
to the relation: C = (ah/e 2 Vk^)(1.05W p ) 1 - 5 /fwj . By 
taking Uu) = 74 meV, we obtain W p = 49(2) meV for 
the LCMO 16 film, and 60(3) meV for the NSMO 16 
film. A larger W p for the NSMO film might be an artifact 
since the residual resistivity p of this film is about 40% 
larger than that of the best single crystal f26|| . The dis- 
crepancy is possibly due to the fact that the interdiffusion 
between the NSMO film and substrate might occur dur- 
ing the high-temperature anneal. From the W p values, 
one can see that for our data T < Wp/ks, which justifies 
the use of Eq. 4. Moreover, the polaron bandwidth is 
greatly reduced compared with the bare bandwidth W Q 
— 4.9 eV. Using g 2 = ln(W /W p ), we can determine g 2 to 
be 4.6(l)/4.4(3) for the LCMO/NSMO 16 film. If hiu 
decreases by 5.7% upon replacing 16 O by 18 O, then our 
calculation shows that the parameter C decreases by 35% 
for LCMO, and by 33% for the NSMO, in good agree- 
ment with the measured values: 35(5)% for LCMO and 
40(7)% for NSMO. 

Furthermore, one can quantitatively explain the iso- 
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tope dependence of E p if small polarons form localized 
bipolarons. In this scenario, <5A = — SW P (see above), 
so 6 A = 0. 057 g 2 W p . From the deduced values for g 2 
and W p above, we calculate that <5A = 12.9 meV for 
LCMO, and 15.0 meV for NSMO. So SE S = 5 A/2 = 
6.5 meV for LCMO, in quantitative agreement with the 
values deduced independently from both resistivity and 
thermoelectric power data. Using 8E p = SE a + SA/2, 
we find 5E p = 13.8 meV for LCMO, and 15.3 meV for 
NSMO. The calculated values are in excellent agreement 
with the observed values: 13.2(3) meV for LCMO and 
14.2(8) meV for NSMO. 

Since the bipolarons in the manganites are immobile, it 
is natural to ask whether bipolarons in the cuprate super- 
conductors, if exist, would be mobile and responsible for 
high-temperature superconductivity [ p7| . The bipolarons 
in cuprates would be mobile if the electron-phonon cou- 
pling in this system were much weaker than in mangan- 
ites. This appears not to be the case since the magnitudes 
of the long-range Frohlich electron-phonon interaction in 
both systems are similar [^lj , and the short-range Jahn- 
Teller (JT) electron-phonon coupling in cuprates is even 
much stronger than that in manganites (in La2Cu04, the 
Jahn- Teller stabilization energy Ejt — 1.2 eV for the Q3- 
type mode @], while E JT ~ 0.5 eV in LaMn0 3 0). The 
bipolarons can become mobile only if the Q3-type "anti- 
JT" mode in cuprates is not active and does not lead 
to the formation of anti-JT polarons possibly due to the 
too large Ejt- Experimentally, the dynamic JT distor- 
tions (or anti-JT polarons) have been observed in doped 
manganites |i~l|-|l4|] , while the dynamic Q3-type JT dis- 
tortions have so far not been observed in doped cuprates. 
Therefore, our present result alone appears not to give a 
positive support to the bipolaronic superconductivity in 
cuprates unless both experiment and theory can clearly 
show that anti-JT polarons (Qa-type) cannot be formed 
in this system. 

In summary, we have observed large oxygen isotope 
effects on the intrinsic resistivity in high-quality epitax- 
ial thin films of Lao.7sCao.25Mn03 and Ndo.ySro.sMnOa, 
and on the thermoelectric power in the ceramic samples 
of Lao.75Cao.25Mn03. The data can be quantitatively ex- 
plained by a scenario [Q where the small polarons form 
localized bound pairs (bipolarons) in the paramagnetic 
state. The coexistence of small polarons and bipolarons 
in the PM state may lead to a dynamic phase separation 
into the insulating antiferromagnetically coupled region 
where the bipolarons reside, and into the ferromagneti- 
cally coupled region where the polarons sit. This simple 
picture can naturally explain the observation of the fer- 
romagnetic clusters in the PM state |p9| . Although we 
cannot completely rule out the possibility that the other 
models might also be able to explain the present isotope 
effects, we would like to point out that the agreement 
between theory and experiment should be quantitative. 
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